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1. Introduction 

The origin of sea quarks of the nucleons remains a subject of intense interest in hadron physics. 
The possible existence of a significant uudcc five-quark Fock component in the proton was pro- 
posed some time ago by Brodsky, Hoyer, Peterson, and Sakai (BHPS) [JlJ to explain the unexpect- 
edly large production rates of charmed hadrons at large forward xp region. The intrinsic charm 
originating from the five-quark Fock state is to be distinguished from the "extrinsic" charm pro- 
duced in the splitting of gluons into cc pairs, which is well described by QCD. The extrinsic charm 
has a "sea-like" characteristics with large magnitude only at the small x region. In contrast, the in- 
trinsic charm is "valence-like" with a distribution peaking at larger x. The presence of the intrinsic 
charm component can lead to a sizable charm production at the forward rapidity (xp) region. 

The CTEQ collaboration [Q] has examined all relevant hard-scattering data and concluded that 
the data are consistent with a broad range of the intrinsic charm magnitude, ranging from null to 2- 
3 times larger than the estimate by the BHPS model. This suggests that more precise experimental 
measurements are needed for determining the magnitude of the intrinsic charm component. 

In an attempt to further study the role of five-quark Fock states for intrinsic quark distributions 
in the nucleons, we have extended the BHPS model to the light quark sector and compared the 
predictions with the experimental data. The BHPS model predicts the probability for the uudQQ 
five-quark Fock state to be approximately proportional to l/ntg, where mg is the mass of the 
quark Q [jl]]. Therefore, the light five-quark states uuduii, uuddd and uudss are expected to have 
significantly larger probabilities than the uudcc state. This suggests that the light quark sector could 
potentially provide more clear evidence for the presence of the five-quark Fock states, allowing 
predictions of the BHPS model, such as the shape of the antiquark x distributions originating from 
the five-quark configuration, to be tested. 

To search for evidence of the intrinsic five-quark Fock states, it is essential to separate the 
contributions of the intrinsic quark and the extrinsic one. Fortunately, there exist some experimental 
observables which are free from the contributions of the extrinsic quarks. As discussed below, 
the d — u and the u + d — s — s are examples of quantities independent of the contributions from 
extrinsic quarks. The x distribution of d — u has been measured in Drell-Yan experiments [^|, |]]. A 
recent measurement of s + s in a semi-inclusive deep-inelastic scattering (DIS) experiment [||] also 
allowed the determination of the x distribution of u + d — s — s. In this paper, we compare these 
data with the calculations based on the intrinsic five-quark Fock states. The qualitative agreement 
between the data and the calculations provides evidence for the existence of the intrinsic light- 
quark sea in the nucleons We also show how the probabilities of various fine-quark states can 
be determined [0]. 

2. Momentum distributions of five-quark states 

For a \uudQQ) proton Fock state, the probability for quark i to carry a momentum fraction jc, 
is given m the BHPS model [0] as 

5 5 m 2 

P(x h ...,x 5 ) =N 5 S(l - J»[mJ - £ ^r 2 , (2.1) 
i=i i=i Xi 
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Figure 1: The x distributions of the intrinsic Q in the uudQQ configuration of the proton from the BHPS 
model. The solid curve is obtained using the expression in Eq. 2.3 for c. The other three curves, correspond- 
ing to c, s, and din the five-quark configurations, are obtained by solving Eq. 2.1 numerically. The same 
probability, ^5 = 0.01, is used for the three different five-quark states. 



where the delta function ensures momentum conservation. N5 is the normalization factor, and nij 



is the mass of quark i. In the limit of >> m pi m\ .2,3, where m p is the proton mass, Eq. 2.1 
becomes 



2 2 

XaXc 



P(x u ...,x 5 )=N 5 !% 2 g(l-I», (2-2) 



where N5 = Ns/rrQ 5 . Eq. 2.2 can be readily integrated over x\, x%, X3 and X4, and the heavy-quark 
x distribution [|lj] is: 

P{xs) = ^V 5 *i[^(l -xs)(l + 10x 5 +4) -2x 5 (l +x 5 )ln(l/x 5 )]. (2.3) 



One can integrate Eq. £3] over X5 and obtain the result = ./V5/36OO, where 2? c ^ is the prob- 
ability for the \uudcc) five-quark Fock state. An estimate of the magnitude of 3?™ was given by 
Brodsky et al. [jl[] as w 0.01, based on diffractive production of A c . This value is consistent with a 
bag-model estimate 

The solid curve in Fig. [I] shows the x distribution for the charm quark (_P(xs)) using Eq. 2.3, 



assuming = 0.01. Since this analytical expression was obtained for the limiting case of in- 
finite charm-quark mass, it is of interest to compare this result with calculations without such an 
assumption. To this end, we have developed the algorithm to calculate the quark distributions us- 



ing Eq. gj] with Monte-Carlo techniques. The five-quark configuration of {xi,...,X5} satisfying 



the constraint of Eq. |2j] is randomly sampled. The probability distribution P(xj) can be obtained 
numerically with an accumulation of sufficient statistics. We first verified that the Monte-Carlo 
calculations in the limit of very heavy charm quarks reproduce the analytical result for P(xs) in 



Eq. 2.3. We then calculated P(xs) using m u = m c j = 0.3 GeV, m c = 1.5 GeV, and m p = 0.938 GeV, 
and the result is shown as the dashed curve in Fig. [[]. The similarity between the solid and dashed 
curves shows that the assumption adopted for deriving Eq. |23| is adequate. It is important to note 



that the Monte-Carlo technique allows us to calculate the quark x distributions for other five-quark 
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Figure 2: Comparison of the d(x) — u(x) data with the calculations based on the BHPS model. The dashed 
curve corresponds to the calculation using Eq. 2.1, and the solid and dotted curves are obtained by evolving 
the BHPS result to Q 2 = 54.0 GeV 2 using fi = 0.5 GeV and fi = 0.3 GeV, respectively. 

configurations when Q is the lighter u, d, or s quark, for which one could no longer assume a large 
mass. 

As shown in Fig. [j], we have calculated the x distributions of the s and d quarks in the BHPS 
model for the \uudss) and \uuddd) configurations, respectively, using Eq. [O. The mass of the 
strange quark is chosen as 0.5 GeV. In order to focus on the different shapes of the x distributions, 
the same value of is assumed for various five-quark states. Figure |l] shows that the x dis- 
tributions of the intrinsic Q shift progressively to lower x as the mass of the quark Q decreases. 
The challenge is to identify proper experimental observables which allow a clear separation of the 
intrinsic light quark component from the extrinsic QCD component. As we discuss next, the quan- 
tities d(x) — il(x), s(x) +s(x), and ii(x) +d{x) — s(x) — s(x) are suitable for studying the intrinsic 
light-quark components of the proton. 

3. Extraction of various five-quark components 

To compare the experimental data with the prediction based on the intrinsic five-quark Fock 
state, it is necessary to separate the contributions of the intrinsic sea quark and the extrinsic one. 
The d(x) — u{x) is an example of quantities which are free from the contributions of the extrinsic 
sea quarks, since the perturbative g — > QQ processes will generate uu and dd pairs with equal 
probabilities and have no contribution to this quantity. The d(x) — u(x) data from the Fermilab 
E866 Drell-Yan experiment at the Q 2 scale of 54 GeV 2 [|[] are shown in Fig. ||. 

In the BHPS model, the u and d are predicted to have the same x-dependence if m u = m^. 
However, the probabilities of the \uuddd) and \uuduu) configurations, and ^l", are not 
known from the BHPS model, and remain to be determined from the experiments. Non-perturbative 
effects such as Pauli-blocking [[J could lead to different probabilities for the \uuddd) and \uuduu) 
configurations. Nevertheless the shape of the d(x) — u(x) distribution shall be identical to those of 
d(x) and u{x) in the BHPS model. Moreover, the normalization of d{x) — u{x) is known from the 
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measurement of Fermilab E866 Drell-Yan experiment [Q] as 

[\d(x)-u(x))dx = ^f-&>l a = 0.m±0m2. (3.1) 
Jo 

Figure^] shows the calculation of the d(x) — il(x) distribution (dashed curve) from the BHPS model, 
together with the data. The x-dependence of the d(x) — u(x) data is not in good agreement with the 
calculation. It is important to note that the d(x) — u(x) data in Fig. |2] were obtained at a rather large 
Q 2 of 54 GeV 2 [Q]. In contrast, the relevant scale, /I 2 , for the five-quark Fock states is expected 
to be much lower, around the confinement scale. This suggests that the apparent discrepancy 
between the data and the BHPS model calculation in Fig. ||] could be partially due to the scale 
dependence of d(x) — u(x). We adopt the value of = 0.5 GeV, which was chosen by Gliick, 
Reya, and Vogt [[!(]] in their attempt to generate gluon and quark distributions in the so-called 
"dynamical approach" starting with only valence-like distributions at the initial \i 2 scale and relying 
on evolution to generate the distributions at higher Q 2 . We have evolved the predicted d(x) — u{x) 
distribution from <2q = ]X 2 = 0.25 GeV 2 to Q 2 = 54 GeV 2 . Since d (x) — u(x) is a flavor non-singlet 
partem distribution, its evolution from <2o to Q only depends on the values of d(x) — u(x) at Qo, 
and is independent of any other parton distributions. The solid curve in Fig. || corresponds to 
d(x) — u{x) from the BHPS model evolved to Q 2 = 54 GeV 2 . Significantly improved agreement 
with the data is now obtained. This shows that the (Revolution should be properly taken into 
account. It is interesting to note that an excellent fit to the data can be obtained if /I = 0.3 GeV is 
chosen (dotted curve in Fig. |J) rather than the more conventional value of /I = 0.5 GeV. We have 



also found good agreement between the HERMES d(x) — u(x) data at Q 2 = 2.3 GeV 2 [11] with 
calculation using the BHPS model. 

We now consider the extraction of the \uudss) five-quark component from existing data. The 
HERMES collaboration reported the determination of x(s(x) + s(x)) over the range of 0.02 < x < 
0.5 at Q 2 = 2.5 GeV 2 from the measurement of charged kaon production in semi-inclusive DIS 
reaction The HERMES data, shown in Fig. [| exhibits an intriguing feature. A rapid fall-off 
of the strange sea is observed as x increases up to x ~ 0. 1, above which the data become relatively 
independent of x. The data suggest the presence of two different components of the strange sea, 
one of which dominates at small x (x < 0.1) and the other at larger x (x > 0.1). This feature is 
consistent with the expectation that the strange-quark sea consists of both the intrinsic and the 
extrinsic components having dominant contributions at large and small x regions, respectively. In 
Fig. H] we compare the data with calculations using the BHPS model with m s = 0.5 GeV. The solid 
and dashed curves are results of the BHPS model calculations evolved to Q 2 = 2.5 GeV 2 using 
/i = 0.5 GeV and jJ. = 0.3 GeV, respectively. The normalizations are obtained by fitting only data 
with x > 0.1 (solid circles in Fig. |f), following the assumption that the extrinsic sea has negligible 
contribution relative to the intrinsic sea in the valence region. Figure [| shows that the fits to the 
data are quite adequate, allowing the extraction of the probability of the \uudss) state as 

&f = 0.024 (/i = 0.5 GeV); &f = 0.029 O = 0.3 GeV). (3.2) 

We consider next the quantity u{x) +d(x) —s(x) —s(x). Combining the HERMES data on 
x(s(x) +s(x)) with the x(d(x) +u(x)) distributions determined by the CTEQ group (CTEQ6.6) [|I"2|], 
the quantity x(u(x) +d(x) — s(x) — s(x)) can be obtained and is shown in Fig. ||. This approach for 
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Figure 3: Comparison of the HERMES x(s(x) +s(x)) data with the calculations based on the BHPS model. 
The solid and dashed curves are obtained by evolving the BHPS result to Q 2 = 2.5 GeV 2 using jU = 0.5 GeV 
and jU = 0.3 GeV, respectively. The normalizations of the calculations are adjusted to fit the data atx > 0.1. 

determining x(u(x) + d{x) — s(x) — s(x)) is identical to that used by Chen, Cao, and Signal in their 



study of strange quark sea in the meson-cloud model Q13[] - 

An important property of U + d — s — s is that the contribution from the extrinsic sea vanishes, 
just like the case for d — u. Therefore, this quantity is only sensitive to the intrinsic sea and can be 
compared with the calculation of the intrinsic sea in the BHPS model. We have 

u(x) + d{x) - s(x) - s(x) = P ua {x u ) + P dd \xj) - 2P s '(xs) . (3.3) 

where P@Q(xq) refers to the .^-distribution of Q in the \uudQQ) state. We can now compare the 
x(ii(x) + d(x) — s(x) — s(x)) data with the calculation using the BHPS model. Since u + d — s — s is 
a flavor non-singlet quantity, we can readily evolve the BHPS prediction to Q 2 = 2.5 GeV 2 using 
/I = 0.5 GeV and the result is shown as the solid curve in Fig. ||. It is interesting to note that a 
better fit to the data can again be obtained with /I = 0.3 GeV, shown as the dashed curve in Fig. ^. 

From the comparison between the data and the BHPS calculations shown in Figs. |2||4], we can 
determine the probabilities for the \uuduu), \uuddd), and \uudss) configurations as follows: 

Sf = 0.122; &>f = 0.240; &f = 0.024 (ju = 0.5 GeV) (3.4) 

or 

@>f = 0.162; 0>f = 0.280; &f = 0.029 (ji = 0.3 GeV) (3.5) 

depending on the value of the initial scale /I. It is remarkable that the d{x) — u(x), the s(x) + s(x), 
and the d(x) + u(x) — s(x) — s(x) data not only allow us to check the predicted ^-dependence of the 
five-quark Fock states, but also provide a determination of the probabilities for these states. 
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Figure 4: Comparison of the x(d(x) + ii{x) — s{x) — s(x)) data with the calculations based on the BHPS 
model. The values of x(s(x) +s(x)) are from the HERMES experiment [|5j], and those of x(d(x) +u(x)) 
are obtained from the PDF set CTEQ6.6 [12|. The solid and dashed curves are obtained by evolving the 
BHPS result to Q 2 = 2.5 GeV 2 using /i = 0.5 GeV and fx = 0.3 GeV, respectively. The normalization of the 
calculations are adjusted to fit the data. 



Equation ^4 shows that the combined probability for proton to be in the \uudQQ) states is 
around 40%. It is worth noting that an earlier analysis of the d — u data in the meson cloud model 
concluded that proton has ~60% probability to be in the three-quark bare-nucleon state [14], in 
qualitative agreement with the finding of this study. A significant outcome of the present work is 
the extraction of the \uudss) component, which is related to the kaon-hyperon states in the meson 
cloud model. It is also worth noting that the \uudQQ) states have the same contribution to the 
proton's magnetic moment as the \uud) three-quark state, since Q and Q in the \uudQQ) states have 
no net magnetic moment. Therefore, the good description of the nucleon's magnetic moment by the 
constituent quark model is preserved even with the inclusion of a sizable five-quark components. 

We note that the probability for the \uudss) state is smaller than those of the \uuduu) and the 
\uuddd) states. This is consistent with the expectation that the probability for the \uudQQ) five- 
quark state is roughly proportional to 1 /nig [jl|, |i5fl . One can now estimate that the probability for 
the intrinsic charm from the \uudcc) Fock state, to be roughly 0.1 &f ~ 0.003. This shows 
that the probability of intrinsic charm could be smaller than the earlier expectation Moreover, 
the <2 2 -evolution would shift the intrinsic-charm distribution to smaller x, suggesting that the most 
promising region to search for evidence of intrinsic charm could be at the somewhat lower x region 
(0. 1 < x < 0.4), rather than the largest x region explored by previous experiments. 



4. Conclusion 



In conclusion, we have generalized the existing BHPS model to the light-quark sector and 
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compared the calculation with the d — u, s + s, and u + d — s — s data. The qualitative agreement 
between the data and the calculations provides strong support for the existence of the intrinsic u, d 
and s quark sea and the adequacy of the BHPS model. This analysis also led to the determination 
of the probabilities for the five-quark Fock states for the proton involving light quarks only. This 
result could guide future experimental searches for the intrinsic c quark sea or even the intrinsic b 
quark sea [16], which could be relevant for the production of Higgs boson at LHC energies [|l7|]. 



This analysis could also be readily extended to the hyperon and meson sectors. The connections 
between the BHPS model, the meson-cloud model [18], the multi-quark models [1^, 20], and the 



lattice QCD calculations [21] should also be investigated. 
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